Accurate computation of Galerkin double surface integrals in the 3-D boundary element method

نویسندگان

  • Ross Adelman
  • Nail A. Gumerov
  • Ramani Duraiswami
چکیده

Many boundary element integral equation kernels are based on the Green’s functions of the Laplace and Helmholtz equations in three dimensions. These include, for example, the Laplace, Helmholtz, elasticity, Stokes, and Maxwell’s equations. Integral equation formulations lead to more compact, but dense linear systems. These dense systems are often solved iteratively via Krylov subspace methods, which may be accelerated via the fast multipole method. There are advantages to Galerkin formulations for such integral equations, as they treat problems associated with kernel singularity, and lead to symmetric and better conditioned matrices. However, the Galerkin method requires each entry in the system matrix to be created via the computation of a double surface integral over one or more pairs of triangles. There are a number of semi-analytical methods to treat these integrals, which all have some issues, and are discussed in this paper. We present novel methods to compute all the integrals that arise in Galerkin formulations involving kernels based on the Laplace and Helmholtz Green’s functions to any specified accuracy. Integrals involving completely geometrically separated triangles are non-singular and are computed using a technique based on spherical harmonics and multipole expansions and translations, which results in the integration of polynomial functions over the triangles. Integrals involving cases where the triangles have common vertices, edges, or are coincident are treated via scaling and symmetry arguments, combined with automatic recursive geometric decomposition of the integrals. Example results are presented, and the developed software is available as open source.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Evaluation of singular and hypersingular Galerkin integrals: direct limits and symbolic computation

Algorithms are presented for evaluating singular and hypersingular boundary integrals arising from a Galerkin approximation in two dimensions. The integrals involving derivatives of the Green's function are de ned as limits from the interior, allowing a simple and direct treatment of these terms. An e cient scheme is obtained by using a combined analytical and numerical approach, the analytic f...

متن کامل

استفاده از روش غنی شده بدون شبکه گلرکین در تعیین پارامتر های شکست صفحات FGM

Stress-intensity factors (SIFs) are the most important parameters in fracture mechanics analysis of structures. These parameters are evaluated for a stationary crack in functionally graded plates of arbitrary geometry using a novel Galerkin based mesh-free method. The method involves an element-free Galerkin method (EFGM), where the material properties are smooth functions of spatial coordinate...

متن کامل

A high precision algorithm for axisymmetric flow

We present a new algorithm for highly accurate computation of axisymmetric potential flow. The principal feature of the algorithm is the use of orthogonal curvilinear coordinates. These coordinates are used to write down the equations and to specify quadrilateral elements following the boundary. In particular, boundary conditions for the Stokes’ stream-function are satisfied exactly. The veloci...

متن کامل

Sinc-Galerkin method for solving a class of nonlinear two-point boundary value problems

In this article, we develop the Sinc-Galerkin method based on double exponential transformation for solving a class of weakly singular nonlinear two-point boundary value problems with nonhomogeneous boundary conditions. Also several examples are solved to show the accuracy efficiency of the presented method. We compare the obtained numerical results with results of the other existing methods in...

متن کامل

Cubature techniques for 3 - D Galerkin BEM

We present cubature methods for the approximation of surface integrals arising from Galerkin discretizations of 3-D boundary integral equations. This numerical integrator is fully implicit in the sense that the form of the kernel function, the surface parametrization, the trial and test space, and the order of the singularity of the kernel functions are not used explicitly. Different kernels ca...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2015